Millikelvin Spatial Thermometry of Trapped Ions 
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We demonstrate millikelvin thermometry of laser cooled trapped ions with high-resolution imaging. This 
equilibrium approach is independent of the cooling dynamics and has lower systematic error than Doppler 
thermometry, with ±5 mK accuracy and =bl mK precision. We used it to observe highly anisotropic dynamics 
of a single ion, finding temperatures of < 60 mK and > 15 K simultaneously along different directions. This 
thermometry technique can offer new insights into quantum systems sympathetically cooled by ions, including 
atoms, molecules, nanomechanical oscillators, and electric circuits. 

PACS numbers: 07.20.Dt measurement of temperature, 37.10.Ty ion traps, 37.10.Jk optical cooling and trapping of atoms 



Laser-cooled trapped ions are a nearly ideal system for 
investigations of quantum physics. The internal ion states 
are strongly decoupled from the surrounding environment 
and can exhibit coherence times of many seconds. In ul- 
trahigh vacuum the trapped ion motion is strongly coupled 
through the Coulomb force but otherwise exhibits good im- 
munity to external perturbations. Precision manipulation of 
ions at the quantum level is readily achieved through the 
use of laser and electromagnetic fields. These properties 
have made laser-cooled trapped ions a preferred platform for 
implementing experiments in quantum information process- 
ing (QIP) |1-4|, and precision metrology (51 IS). The strong 
Coulomb coupling makes laser-cooled trapped ions attractive 
for sympathetic cooling at millikelvin temperatures in investi- 
gations of fundamental physics |7 |, the dynamics of complex 
molecular 111 and biomolecular |9| ions, nano-mechanical 
oscillators fTO", TT|, resonant electric circuits [11211 . and Bose- 
Einstein condensatesi 13 |. 

Thermometry is a key diagnostic in all of these exper- 
iments. Effective Doppler cooling is a precondition for 
QIP and precision measurement applications. The devel- 
opment of new ion trapping architectures for QIP requires 
millikelvin thermometry to quantify high anomalous heating 
rates |[T4l [T5]| . Motional decoherence is a limiting factor in 
recent experiments including the long range mechanical cou- 
pling of ionsi 16 , 17 1 and the realization of a phonon laser ifTSl . 
In sympathetic cooling, thermometry of the cooling ions is the 
primary diagnostic for the temperature of the cooled species. 
In the presence of a controlled heating process, thermometry 
of the cooling ions can be used to measure the cooling power 
of the sympathetic cooling process. Sympathetic cooling is a 
necessary precondition for mixed species experiments such as 
the Al+/Be+ ion clock |5 1 or molecular ion experiments! 7- 
|9l. These applications of thermometry are especially useful 
for proposed hybrid systems in which trapped ions sympa- 
thetically cool and interact with nanomechanical oscillators 
or electric circuits. 

Existing trapped ion thermometry techniques in the mil- 
likelvin regime are spectroscopic and rely on varying the 
cooling laser frequency to obtain a measurement of Doppler 
broadening due to ion motion |[T9l[2Ql . This approach pushes 
the system out of thermal equilibrium, with the variation in 




FIG. 1: Schematic of the experimental apparatus. A single Yb^ ion 
is confined at the node of an RF electric quadrupole (dashed lines) 
formed between two needle electrodes. The ion is laser cooled using 
a A = 369.5 nm laser beam (arrow) and the resulting fluorescence 
is collimated using an in- vacuum phase Fresnel lens. This light is 
subsequently imaged onto a cooled CCD camera (not shown). 

laser frequency affecting the cooling rate and subsequently 
changing the temperature over the course of the measurement. 
Recoil heating of the ion increases its temperature and re- 
duces the peak fluorescence when the laser detuning is less 
than half the linewidth from resonance. At laser detunings 
larger than half the linewidth the cooling rate is reduced caus- 
ing an increase in the temperature and the associated Doppler 
linewidth. Both of these effects produce a systematic error 
that leads to an overestimation of the temperature. Time re- 
solved techniques to quantify the effects of non-equilibrium 
laser cooling dynamics on the fluorescence signal have been 
demonstrated |[T5l [2T1l but depend on numerous simplifying as- 
sumptions or extensive simulation! 22 1. The Doppler approach 
is also limited to measurement of the ion motion in a sin- 
gle dimension fixed along the axis of the laser beam. The 
spectroscopic linewidth of the transition limits the minimum 
detectable temperature by the Doppler technique. These defi- 
ciencies are less of an issue in traps with low anomalous heat- 
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ing rates fT4l where sub-millikelvin resolved sideband ther- 
mometry is feasible. 

We demonstrate an alternative approach to trapped ion ther- 
mometry which probes the spatial distribution rather than the 
velocity distribution. A key enabling technology of this tech- 
nique is our recent demonstration of high resolution imaging 
of trapped ions with spot sizes as small as 373 nm ra- 
dius (440 nm FWHM) |23|. Such high resolution enables 
us to measure temperatures of a few millikelvin by spatial 
thermometry and investigate phenomena inaccessible through 
spectroscopic thermometry. For a particle in a harmonic po- 
tential the spatial extent is dependent on the temperature and 
the trapping frequency. Imaging an ion thus provides a mea- 
sure of its temperature in two dimensions and under steady- 
state conditions. Other trapped ion systems typically have 
imaging resolutions of several /im and are thus limited to a 
spatial temperature sensitivity of a few Kelvin. Our accuracy 
of ±5mK is limited by uncertainty in the imaging resolution 
while the precision of ±lmK is limited by image signal to 
noise. 

The apparatus is similar to that described previously |[23l. 
174y]3+ ions were generated by isotope- selective pho- 
toionisation of a Yb beam and were trapped in an RF 
electric quadrupole field formed by applying a potential 
Vocos{^}RFt), Qrf/^tt = 20 MHz between two needles 
separated by 300 /im. We obtained ion motional frequen- 
cies between 500 kHz and 1500 kHz. Residual electric fields 
at the trap center were canceled by applying DC voltages to 
compensation electrodes. The ions were laser cooled using 
A = 369.5 nm light resonant with the S1/2 to P1/2 transi- 
tion. The frequency of the cooling laser was stabilized using 
dichroic atomic vapor laser locking (DAVLL) |24|. Detuning 
the laser frequency changed the DAVLL offset which in turn 
could be calibrated using a Fabry-Perot interferometer with a 
known free spectral range. The ion fluorescence was collected 
using a phase Fresnel lens (PFL) mounted in the ultra-high 
vacuum chamber, close to the trapping region |25|. The PFL 
has an numerical aperture (NA) of 0.64 which covered 12% of 
the solid angle and collected 4.5% of the total emission. The 
collimated fluorescence was imaged with a magnification of 
596 ± 22 on to a cooled CCD camera. 

We initially made a qualitative investigation of the spatial 
thermometry technique by controllably heating the ion. To 
heat the ion, we applied voltage with a white noise spectrum 
and a bandwidth encompassing all motional frequencies of the 
ion to one of the compensation needles. This allowed us to in- 
vestigate different temperature ranges in a controlled fashion. 
Figure [2] shows ion images for three different levels of heat- 
ing; (a) no external heating, (b) low external heating and (c) 
high external heating (3.3x the noise voltage of b). There is 
a clear increase in the ion spot size with the increase in ion 
temperature. 

From these high resolution images, we can quantitatively 
calculate the temperature of the ion. We fit a two dimensional 
Gaussian to the ion image in order to quantify the ion im- 
age radii Wx^y for horizontal (x) and vertical (y) axes. 
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FIG. 2: Images of a trapped ion for three different external heating 
rates. Temperatures were calculated using equation [T] from the hor- 
izontal (laser cooling direction) ion spot radii. In the case of no (a), 
low (b) and high (c) heating rates we calculate temperatures of 5 =b 5 
mK, 30 =b 9 mK and 162 zb 14 mK. Errors were dominated by sys- 
tematic uncertainty. 

The vertical axis of the image is parallel to one of the trap 
axes. The motional frequencies along the other two trap axes 
are degenerate. Assuming a thermal distribution 1261 , the ion 
temperature along the x axis is given by 

Tx = (1) 

Kb 

where m is the mass and /c^ is Boltzmann's constant and (x^) 
is the variance of the ion spatial distribution along the x axis. 
Owing to the degeneracy of the other two trap axes, a similar 
equation holds for y. 

In the low temperature limit, the finite resolution of the 
imaging system dominates the observed ion image radius. Us- 
ing our knowledge of the imaging system, we can estimate the 
variance of the ion spatial distribution and hence the tempera- 
ture. To extract the variance of the ion spatial distribution we 
deconvolved the measured ion image radius with the imaging 
resolution (wi) according to {x ) = '^^ . Uncertainty 
in the imaging resolution leads to a systematic uncertainty in 
the variance of the ion spatial distribution and thus the mea- 
sured temperature. We evaluate the uncertainty by consider- 
ing upper and lower bounds on our imaging resolution. An 
ideal imaging system has a minimum 1/e^ radius resolution 
of 0A3X/NA. In our system this gives a diffraction limited 
lower bound on the imaging resolution of 249 nm. The small- 
est spot radius observed in the horizontal axis was 373 nm, 
giving an experimental resolution limited upper bound. De- 
convolving the radii of the ion image shown in Fig [2] a) with 
the lower resolution bound gives a temperature of 8 ± 1 mK. 
Here the statistical error of 1 mK arises from the error in es- 
timation of Wx and sets the resolution limit of spatial ther- 
mometry. The spot size of this image is the smallest we have 
observed, so deconvolution with the upper bound of imaging 
resolution gives a temperature estimate of approximately zero. 
The deconvolution procedure significantly improved the accu- 
racy of our technique in the few mK temperature range. This 
systematic uncertainty is the primary limit to the accuracy of 
spatial thermometry. In the following, we state systematic un- 
certainties for imaging temperatures unless indicated other- 
wise. 

We compared the imaging thermometry technique to the es- 
tablished Doppler spectroscopic thermometry technique. The 
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ion fluorescence spectrum under controUed heating was mea- 
sured by sweeping the laser frequency vl over the cooHng 
transition for various levels of applied noise voltage. When 
the laser was detuned to frequencies above resonance, the ion 
was rapidly heated causing a sharp decrease in the fluores- 
cence. The half Voigt profile that is conventionally used to 
fit the spectral response 1 27 1 was not able to discriminate be- 
tween Lorentzian and Gaussian linewidth contributions. In- 
stead we fit the data to a Lorentzian with a smoothed step 
function cutoff 
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where 5vl is the laser frequency detuning, Ajr^ is the linewidth 
of the cooling laser and Tt is the FWHM spectral linewidth of 
the ion. To determine the temperature using the spectroscopic 
thermometry technique, the total spectral linewidth (Tt) was 
deconvolved into its two constituent components |28|. The 
Gaussian component (Tg) is related to the temperature of the 
ion by T = (FgA)^ |29 | and the Lorentzian component 
(Fl) is related to the linewidth of the ion at zero temperature. 

A systematic uncertainty in the temperature estimated by 
the spectroscopic thermometry technique arises from unre- 
solved contributions to the Lorentzian component of the spec- 
troscopic linewidth. We calculate the upper bound for the 
spectroscopic temperature by assuming that F^ is the natural 
linewidth, resulting in a maximum contribution from F^. This 
bound assumes that broadening mechanisms such as Zeeman 
splitting, saturation and micromotion have no contribution to 
F L . The lower bound on the spectroscopic temperature is cal- 
culated by assuming F^ is equal to the smallest experimen- 
tally observed linewidth. 

Figure [3] compares the spectroscopic temperature results to 
those from the imaging method. The same three external heat- 
ing rates as in Fig |2] are used and imaging temperatures are 
obtained for several cooling laser detunings. The imaging 
temperatures at each detuning were calculated from the fit- 
ted ion spot sizes in the horizontal (cooling laser) direction. 
In all cases the imaging thermometry is far more precise and 
appears to be more accurate. The temperatures calculated by 
spectroscopic thermometry show large systematic uncertain- 
ties in all cases. The spectroscopic thermometry technique 
overestimates the temperature as expected from considera- 
tions of laser cooling dynamics. 

Spatial thermometry shows dynamics that lead to system- 
atic error in spectral thermometry. At large detunings an 
increase in temperature was observed as the scattering rate 
dropped, reducing the efficiency of Doppler cooling. When 
the cooling laser was close to resonance an increase in the 
temperature occurred due to increased recoil heating. In Fig 
[3] a), with no applied heating, a minimum in the ion imaging 
temperature was observed at ~ 15 MHz corresponding to the 
expected optimum detuning of Ft/2. The systematic uncer- 
tainty in the imaging temperature dominates the measurement 
accuracies for temperatures close to the Doppler cooling limit. 
As external heating is applied and the temperature of the ion 
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FIG. 3: Ion temperature dependence on cooling laser detuning for 
three different external heating rates. Both spectroscopic and imag- 
ing temperatures for (a) no, (b) low, and (c) high external heating 
are shown. The uncertainty in the imaging temperature is dominated 
by the systematic uncertainty in the correction for ion spatial extent 
by the imaging resolution. The spectroscopic temperature la uncer- 
tainty band for each of the three external heating rates is shown as 
the gray band with the upper and lower bound depicted as a dashed 
line. 



increases, the systematic uncertainty becomes negligible. In 
Figure |3](c), imaging temperatures for detunings greater than 
70 MHz were not recoverable due to the configuration of the 
imaging system. 

The spatial thermometry technique allows measurement of 
temperature along two dimensions even when the ion dynam- 
ics are highly anisotropic. Figure]?] shows the effect of highly 
anisotropic laser cooling dynamics for a single ion. The cool- 
ing laser was controllably decoupled from one axis of motion 
by applying a DC bias to the RF needles, which caused ro- 
tation of the trap axes. When the cooling laser is decoupled 
from one axis of motion, the cooling power is reduced while 
there is no reduction in recoil heating along this axis. This ro- 
tation is believed to occur due to the difference in the RF and 
DC grounding. Figure |4] shows the temperature of the ion as 
a function of trap axis rotation demonstrating a large disparity 
between the horizontal and vertical ion temperatures. Ion im- 
ages in Fig [4] show the ion spatial extent at selected rotations. 
A temperature difference of a factor of 1000 was observed 
without any adverse effects on the trapping. This temperature 
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difference was achieved with a maximum trap aspect ratio of 
2.25. The angular width of the feature in Fig|4]b) was much 
larger than the angular width of the beam used to cool the ions 
(0.3°). Spectral thermometry along the vertical axis would re- 
quire the laser direction to have a projection along that axis, 
inducing laser cooling. Hence, the anisotropic dynamics stud- 
ied here can only be observed by spatial thermometry. 
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FIG. 4: Spatial thermometry of a single ion under highly anisotropic 
laser cooling dynamics. Top: ion images for several rotation angles. 
The contrast of the ion images has been adjusted for ease of viewing. 
Bottom: temperature along weak-cooling axis (filled black squares) 
and strong-cooling axis (open red circles). For rotation angles be- 
tween -3.5° and -2.5° (gray shaded area) the size of the ion image is 
larger than the camera chip along the vertical direction. A tempera- 
ture difference of a factor of 1000 was sustained without any adverse 
effects on the trapping. Error bars are smaller than the data points. 

High-resolution imaging of trapped ions was used to realize 
steady state millikelvin thermometry and investigate phenom- 
ena inaccessible through conventional Doppler spectroscopy. 
We achieved an accuracy of ±5 mK with temperature resolu- 
tions of ±1 mK, limited by the imaging system resolution. 
This spatial thermometry technique is independent of laser 
cooling dynamics and has been used to observe the depen- 
dence of ion temperature with laser detuning for three differ- 
ent heating rates. The variation of temperature with cooling 
laser detuning is a prominent source of systematic error in 
Doppler thermometry and plays an important role in estimat- 
ing heating rates [ 22 1 . By rotating the trap axes with respect 
to the cooling laser we investigated the ion dynamics under 
highly anisotropic laser cooling. Temperatures of < 60mK 
and >15K were simultaneously measured along two motional 
axes of a single trapped ion. The latter temperature is eas- 
ily accessible by optical cryostats and demonstrates the util- 
ity of this technique for investigations of sympathetic cooling 
in proposed hybrid trapped ion systems with nanomechani- 
cal oscillators 1101 lUJ or resonant electrical circuitsi 12|. In- 
vestigations of highly localized variations in anomalous heat- 
ing rates observed in micro-fabricated ion traps developed for 



QIPfTTl suggests a strongly anisotropic process that could be 
better resolved using the multi-dimensional aspects of spatial 
thermometry. 
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